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ABSTRACT 

We have investigated the role of photometric variabiUty in causing the apparent age 
spreads observed in the colour-magnitude diagrams of OB associations. We have found 
that the combination of binarity, photometric uncertainty and variability on timescales 
of a few years is not sufficient to explain the observed spread in either of the OB 
associations we have studied. Such effects can account for about half the observed 
spread in the cr Orionis subgroup and about 1/20 of the observed spread in Cep 
0B3b. This rules out variability caused by stellar rotation and rotation of structures 
within inner accretion discs as the source of the majority of the the apparent age 
spreads. We also find that the variability tends to move objects parallel to isochrones 
in V/V — i' CMDs, and thus has little influence on apparent age spreads. We conclude 
that the remaining unexplained spread either reflects a true spread in the ages of the 
PMS objects or arises as a result of longer term variability associated with changes in 
accretion flow. 
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1 INTRODUCTION 

A significant unresolved question in the study of star for- 
mation is how long it takes. Whether we are considering the 
entire process from the state of neutral interstellar hydrogen 
to the Zero Age Main-Sequence (ZAMS), or just the portion 
of the process from the fragmentation of a giant molecular 
cloud (GMC) onwards, there is no consensus as to how long 
star formation takes. There are two competing paradigms 
of star formation currently proposed in the literature, which 
each give rise to very different timescales for star formation. 

1.1 Slow star formation 

IShu ri977t) put forward the model of star formation which 
we will re fer to as slow star formation (SSF), which was re- 
viewed bv lShu et al. ('1987^ . In this model, molecular cloud 
complexes are in dynamical equilibrium, with lifetimes of 
several tens of megayears and are supported against free fall 



collapse by magnetic fields. However, in clumps where the 
density is n ^ 10^ cm~'^, the ionisation fraction can be suffi- 
ciently low that the neutral molecular material may diffuse 
through the magnetic field, removing flux from the core (am- 
bipolar diffusion), leading to higher degrees of central con- 
densation. Eventually the field can no longer support the, 
now prestellar, core against gravitational collapse to form a 
hydrostatic protostar. 

Typically, predicted ambipolar diffusion timescales lie 
in the range 5-10 Myrs. This is consistent with the observed 
lifetimes of pre-stellar cores w ith 10^ < n < 10^ found by 
IWard- Thompson et al. to be ~ 10 yrs, through com- 

parison of the number of star less cores versus protostellar 
cores. However, Ijiiina et al. (1 999) performed a more com- 
prehensive study of pre-stellar cores using ammonia emis- 
sion. They found the ratio of starless to stellar cores to be 
too small to be consistent with that expected if the ambipo- 
lar diffusion timescale governed core lifetimes. This suggests 
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that the timescale of the starless core phase is much less 
than the ambipolar diffusion timescale. 

1.2 Rapid star formation 

More recently, an alternative paradigm for pre-stellar core 
form ation and collapse has been discussed in the literature 
fe.g. iBa llesteros-Paredes et al. 19991 : iHartmann et al. 200lt 
iHartmann 200ll : lElmeereen 20001)" In this model the star 
formation rate is regulated by supersonic turbulence, 
not magnetic fields. Since supersonic turbulence in self- 
gravitating clouds is expected to decay rapidly, the cloud 
support mechanism in this pictur e leads to a short cloud 
lifetime compared to that in SSF dPrinde et al. 200 ill . So, 
rather than being long-lived structures, GMCs are treated 
as transient objects that form through interactions between 
supersonic flows in the interstellar medium. The whole pro- 
cess from formation of the GMCs through to the arrival 
of the pre-main-sequence (PMS) stars on the birthline is 
expected to take around 3 Myrs, with collapse of cores 
to form protostars occurring almost immediately. As such 
we will refer to this model as rapid star formation (RSF). 
RSF has a number of distinct advantages over SSF, both 
observationally and theor etically which are reviewed by 
iMac Low fc Klessen (200411 . 

1.3 Age spreads 

The presence of apparent age spre ads has been observed 
in a number of young clusters (e.g. iHerbst fc Miller 1982 ; 
Sung_et_aLJ.99^ and associations fe.g. |Pozz^^^d"^003 ; 
Dolan fc Mathieu 20011 ) and their sizes have been pre- 
sented a s evidence in favour of both SSF and RSF. 
IPalla fc"stahler (2000) investigated apparent age spreads in 
a number of young star forming regions and found evidence 
of accelerating star formation, which they use to argue in 
favour of SSF. A further study by the same authors into the 
age spread in the Taurus- Auriga region also fo und ev idence 
of accelerating star formation (Palla & Stahlgr, 2002h . The 
evidence of accelerating star formation in the regions studied 
by these authors came in the form of distributions of stel- 
lar ages that were, generally, strongly pea ked at 1 - 2 Myrs . 
with few PMS stars older than this. As iHartmarm (20031) 
pointed out, this surely implies there is something special 
about the last 1-2 Myrs, if such widely separated regions 
have formed the majority of their stars at the same time, 
whilst their overall lifetimes are ~ 10 Myrs. Such an uncom- 
fortable, special, state of affairs is not required if the one 
accepts the RSF paradigm. This is because the strong peak- 
ing of age distributions at 1-2 Myrs follows naturally if this i s 
the timescale for cloud and star formation (|Hartmann 2003() . 

ElmcErcen (2000) found that the ages spreads seen in 
a number of OB associations and clusters are comparable 
to the inferred crossing time for the parent cloud. The con- 
clusion that is drawn from this is that the age spread is in- 
dicative of the timescale for star formation, and that this is 
comparable to the crossing time, as expected for RSF. How- 
ever, his data reveal that the age spreads also scale with 
the mean ages for the groups he uses. This scaling of appar- 
ent age spreads with mean age for clusters and associations 
suggests that the age spreads originate from a photometric 



scatter of given magnitude. The size of this spread is much 
larger than any photometric uncertainties. Since PMS ob- 
jects move more slowly through colour-magnitude space at 
older ages, a given photometric scatter will naturally imply 
larger age spread as the mean age rises. The correlation of 
spread with inferred crossing time may also arise as a re- 
sult of the fact that older open clusters and associations are 
larger and thus a longer crossing time is inferred for the 

par ent cloud. 

IHartma nn (20oill 20031) disputed th e results and con- 
clusions of IPalla fc Stabler (2000l 2002h . arguing that the 
actual age spreads are much smaller than those observed, 
but accepted that the presence of such large age spreads 
(10^ yrs) would be a problem for the RSF paradigm. 
Additionally, it is not clear that the interpretation of 
the age spreads being indicative of timescale is correct. 
iTassis & Mouschov ias (2004) pointed out that interpreta- 
tion of age spreads in this manner already assumes a core 
formation timescale that is essentially instantaneous with 
respect to the lifetime of the molecular cloud, r„ic. In real- 
ity, an age spread could only tell us the difference between 
Tmc and the timescale for forming a body that will survive 
the destruction of the molecular cloud, the core formation 
timescale, Tc/. So: Tspread = Tmc — Tcf- As such, interpre- 
tation of any spread in ages is dependent on knowledge of 
either r^/ or 

Incidentall y, this point is also relev ant to the argument 
put forward bv IHartmann et al. (20011) that, since very few 
young clusters or associations older than 3 Myrs are associ- 
ated with their parent cloud, star formation must occur on 
a timescale of about 3 Myrs. This also implicitly assumes 
that Tcf is short. The same observational evidence could in- 
dicate that star formation is slow, but that the cloud is dis- 
rupted quickly after the first PMS stars arrive at the birth- 
line. However, as pointed out by Hartmann ct al. (200]J), the 
vast majority of local molecular cloud complexes show evi- 
dence of star formation. This must imply that r^f should be 
short, or the question must surely be asked as to the where- 
abouts of the clouds which are still mostly in the pre-stellar 
phase of star form ation. This is reflected in the results of 
Ijiiina et al. (1999h . described earlier, that the ratios of pre- 
stellar to stellar cores found in molecular clouds are far below 
the 3 to 30:1 range required by SSF. 

It is important to recognise, however, that the reality 
of apparent age spreads has not been well established. The 
principal evidence for spreads of ages actually comes from 
an observed spread of PMS stars in colour-magnitude (C-M) 
space. A single age would be expected to give rise to a much 
more narrow distribution of stars about an isochrone. There 
are a number of plausible alternatives that could explain 
such C-M spreads for a population that in reality arrived at 
the birthline simultaneously. For example accretion-driven 
age spread (^out ct al. 1999} could give rise to an appar- 
ent spread of ages for an ensemble o f objects with differing 
accretion histories. iHartmaim (20011) explored a number of 
other factors that might be expected to influence the degree 
of observed age spread. These included variable extinction, 
photometric variability, differences in accretion luminosities 
and the presence of unresolved binaries. Establishing the in- 
fluence of variability and binarity on the width of the PMS 
in OB association-like environments is the aim of this work. 
We have obtained 2 epoch, 2 colour photometry for 2 regions 
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with differing photometric spreads to investigate this effect. 
The 2 regions studied here are within weii known OB as- 
sociations: Cep OB3b and the a Orionis young group (part 
of the Orion OBlb association). We have simulated the de- 
gree of spread introduced by variability on timescales of less 
than 1 year in Cep OB3b and less than 4 years in a Ori. By 
comparing 2 colour catalogues of PMS objects within these 
associations we have estimated their variability, and used 
this to simulate the PMS in C-M space, assuming a single 
isochronal age for the objects. 

The rest of the paper is laid out as follows. In section |5| 
we describe our observations and data reduction. In section 
l^we describe how we have simulated the spreads in each of 
our associations in turn, and give the basic results. These 
results are discussed in section |1] and our conclusions are 
summarised in section |5] 



the residuals to the 6-coefRcient fits were all less than 0.1 
arcsec. 



2 OBSERVATIONS AND DATA REDUCTION 

All observations were carried out using the Wide Field Cam- 
era (WFC) mounted on the Isaac Newton Telescope (INT) 
at the Roque de los Muchachos Observatory, La Palma. Our 
first epoch data set f or a Ori is the same as that presented by 
iKenvon et al. (200^ . It was made up of observations taken 
on the nights of 27-30 September 1999 of 5 fields of view 
(FoV) in the Harris R and Sloan i' filters. We carried out 
observations of 4 FoVs coincident with Kenyon et al's sur- 
vey on the night of 7 September 2003, using the same filters. 
The new observations are detailed in Table Q We have ob- 
tained new data for both epochs for our Cep OB3b survey. 
We have observed a single WFC FoV in this region, chosen 
to cover the area with the highest density of PMS objects 
identified by Pozzo ct al. (2003). Observations were taken 
on the nights of 12 September 2003 and 28 September 2004 
using Harris V and Sloan i' filters. Again the observations 
are detailed in Table 

Data obtained in 2003 and 2004 were reduced 
in an identical manner. Flatfields and data frames 
were linearised using the 2003 August coefficients (see 
[http:/ / www.ast.cam.ac.uk/~wfcsur/technical/foibles/index.phpl 
for details), and then bias subtracted using a median bias 
frame specific to each night. We flat fielded the data using 
frames constructed from twilight sky flatfields taken during 
the same observing runs as the data being corrected. The 
i'-band frames were successfully defringed using a fringe 
frame from 2001 September, obtained from the web pages 
of the Cambridge Astronomy Survey Unit (CASU). 

Optimal photometry w as pe rformed using the m ethod 
laid out by iNavlor (19m) and iNavlor et al. (2003). with 
the revisions described in iBurningham et al. (2003) and 
iLittlefair ct al. (2005). We have allowed the profile correc- 
tion to vary with position on each chip, and fitted it with a 
3rd order polynomial in the x-axis and a 5th order polyno- 
mial in the y-axis. We do not apply the 'ill-determined sky' 
fiag for the Cep 0B3b data, and we allow the thresholds for 
its application to vary from field to field for the a Ori data. 
For obvious reasons we did not apply the 'variable' data 
quality flag to reject objects that showed evidence of vari- 
ability. An astrometric solution was obtained through com- 
parison with a 2MASS catalogue for each FoV. The RMS of 



As described by INavlor et al. Hool, overlap regions 
between pointings were used to normalise the catalogues for 
each pointing onto a single system. We use the same method 
to bring observations from 2 epochs for each of our target 
regions to the same system. In the case of the a Ori data, the 
way in which flatfields were normalised had changed between 
the two epochs, so the data from each epoch were treated 
slightly differently. When the 1999 data were reduced, the 
fiatfield for each chip on the WFC was normalised sepa- 
rately. By the time the 2003 data were reduced, the reduc- 
tion software had been changed such that the fiatfields for all 
4 chips in each pointing were normalised together. As such, 
prior to catalogue normalisation, the 1999 exposures were 
combined to produce 1 catalogue for each chip, whilst the 
2003 exposures were combined to produce 1 catalogue for 
each pointing. The overlap regions between the interlocking 
pointings were then used to normalise the catalogues onto a 
single system. 

To verify that there were no major sources of error intro- 
duced at the image processing stage we checked that the dis- 
tribution of positive and negative differences between obser- 
vations were uniform with position on the sky for the a Ori 
data. We found that the trends in the differences between 
observations were correlated with the pixel coordinates on 
each chip from the 1st epoch catalogue. iKenvon et al. (20051) 
found the RMS of differences between overlap stars seen in 
their data to be 0.05 mags and speculated that it was due 
to problems with the flatfield, which is consistent with what 
has been seen here. They state that the camera was suffering 
from light leaks during their observing run, and it is likely 
that this is the origin of the problem. We fitted the trends 
in the differences in x and y for each chip with 1st and 2nd 
order polynomial functions, and corrected the magnitudes 
of the objects in the 1st epoch catalogue accordingly, prior 
to re-normalising the 2 epochs. The RMS of the differences 
in the magnitudes in the interlocking overlap regions after 
normalisation suggested that an additional magnitude inde- 
pend ent uncertainty of a bout 0.01 mags was also present. 
As in iNavl or et al. r2002h . this uncertainty was included in 
the uncertainty estimate for all objects. Following this we 
found that positive and negative differences for the a Ori 
data were spatially uncorrelated. 

Since no overlap regions were present in the Cep 0B3b 
data sets, we could not use the results of the normalisation 
to determine the size of any additional uncertainty, so we 
applied the same value as found for the a Ori data set. The 
distribution of positive and negative differences were found 
to be spatially uncorrelated in the Cep OB3b data set. 

We have not applied a measured photometric calibra- 
tion to our final catalogues as this would have involved trans- 
forming our Harris and Sloan magnitudes and colours into 
the Cousins system. Such a transformation would risk intro- 
ducing spurious correlations between variability and colour, 
whilst not improving our experiment in any way. As such we 
have simply applied mean zero points to our magnitudes and 
colour coefficients to our colours to construct our CMDs. 
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Table 1. Summary of Observations 



Field Name 


a Ori 1 


cr Ori 2 


(7 Ori 3 


(T Ori 4 




Cep OB3b 


RA(J2000) 


05 40 14.2 


05 40 13.1 


05 38 07.7 


05 38 07.4 




22 55 43.3 


Dec(J2000) 


-02 20 18.1 


-02 51 48.0 


-02 20 18.0 


-02 51 51.0 




-1-62 40 13.7 


Filter 


R i' 


R i' 


R i' 


R i' 


V 


i' 


2003-09-07 


600s 300s 


600s 300s 


600s 300s 


600s 300s 






2003-09-12 










30s, 


2s 13s, 2s 


2004-09-28 










30s, 


5s 300s, 30s, 5s 



3 SIMULATION OF PHOTOMETRIC 
SPREADS 

3.1 The (T Ori subgroup 

The o Ori young group has been the subject of a num- 
ber of spectroscopic studies aimed at i dentifying bona fide 
members. In particular two campaigns l)Kenvon et al. 200,4 
[Bur ningham et al. 20 05i) were successful in identifying a 
large number of members present in the catalogue obtained 
from the 1st epoch observations. Importantly, these spec- 
troscopic surveys also demonstrated that photometric selec- 
tion of the PMS objects in this region is not subject to seri- 
ous contamination, and does not miss significant numbers of 
bona fide members. As such we have used the membership 
lists from these studies to define a PMS region in the CMD, 
with minimal contamination, from which we have drawn our 
sample (see Figure 0. Only objects with little doubt as to 
their membership have bee n used to guide the photo metric 
selection. In the case of the iBurningham et al. ^2005^ cata- 
logue this meant objects with greater than 90% membership 
velocity probability. In the case of the Konyon ct al. f2005|) 
catalogue this meant objects which displayed strong Li ab- 
sorption, evidence of low surface gravity and the appropriate 
radial velocity. 

To test if the observed C-M spread is affected by vari- 
ability we have assumed a single age for the young group 
and then simulated the PMS using an estimate of each ob- 
ject's variability, derived from our 2 epoch observations. This 
method implicitly includes the effect of photometric uncer- 
tainty. 

To verify that our PMS selection is indeed more vari- 
able than the background we have compared the RMSs of 
the differences between the two sets of observations for the 
background objects and for the PMS selection. To make such 
a comparison meaningful we have restricted our sample to 
those i' magnitudes where the uncertainties are small. In 
Figure |21 we plot the uncertainty in the R — i' differences 
(thus incorporating uncertainties from both epochs) against 
i' (2nd epoch). As can be seen, the uncertainties rise sharply 
fainter than i' = 18, as the uncertainty in the photon count- 
ing statistics of the sky begins to dominate. Brighter than 
i' = 16, some objects also display higher uncertainties. This 
is because for these points the 1st epoch data were drawn 
mainly from a short exposure. In the range 15 < i' < 18 the 
uncertainties are dominated by the small magnitude inde- 
pendent uncertain ty measured during the normalisation of 
the catalogues fsee lNavlor et al. 20021 and Section|5J . Based 
on this plot we restrict our comparison of the RMSs to the 
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Figure 1. The CMD for the 2nd epoch Ri' cata- 
logue. Spectroscopic m embers of iKenvon et al. f2005^ and 
IBurningham et al. (20051) are shown as asterisks. Objects iden- 
tified as non-members in the same studies are shown as open 
circles. The solid line indicates the location of the cut for our 
PMS selection. 



range 15 < i' < 18. The RMS of the differences between the 
two sets of observations indicate that the objects in the pho- 
tometric PMS selection are significantly more variable than 
those in the background. The RMS for the differences in i' 
are 0.05 for the background and 0.09 for the PMS region. In 
R-i' the RMSs are smaller: 0.02 for the background; 0.04 
for the PMS region. 

We placed the objects on a single empirical isochrone 
by fitting a 3rd order polynomial to the photometrically se- 
lected PMS in the CMD. Each point used for fitting the 
polynomial is weighted according to its uncertainty in R — i' , 
with the points carrying the greatest uncertainty having the 
smallest influence on the fit. A point was then placed on 
this empirical isochrone at the appropriate magnitude for 
each PMS object. One can view this as moving each point 
in colour until it falls on the isochrone. 
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Figure 2. A plot of the uncertainty in A{R — i') against i' for 
the FoV centered on the a Ori subgroup. 



Figure 3. The mean xi the simulated PMS compared to 
the observed PMS plotted against binary fraction, shown as a 
percentage for the PMS objects near cr Ori. 



Having placed each point on the isochrone, we simulated 
the effect of binaries by splitting this empirical isochrone 
into a single star and an equal mass binary sequence. We did 
not assume a binary fraction, but rather simulated spreads 
with a range of fractions. For each binary fraction we ran- 
domly selected that proportion of objects to be offset in 
magnitude by an amount equal to 0.75(1 — f tin), and the rest 
by 0.75 — 0.75(1 — fbin)- This gave rise to a difference of 0.75 
mag between the two sequences, as would be expected from 
equal mass binaries. We do not attempt to simulate q < 1 
binaries as these would fall within the envelope bounded by 
our two sequences, and thus would not have a significant 
impact the results of this investigation. 

To simulate the effect of variability we next moved each 
point by an amount in both magnitude and colour equal to 
{A/y/2), where A is the observed difference between the two 
observations for that point. Characterising the variability in 
this manner for each object in turn has distinct advantages 
over parameterising the scatter of the whole sample, as any 
correlations between Ai' and A{R — i') have been included 
without any assumptions as to the source of the variability. 
Furthermore, any correlation of variability with magnitude 
has also been included, as has the influence of photometric 
uncertainty. 

A further source of scatter that might influence the ob- 
served spread is differential reddening. We have neglected 
the effects of differential reddening in this simulation since 
the extinction towards a Ori is low , with a colour excess of 
just E(B -V) = 0.05 ijLee 196ft) . To select which binary 
fraction gave the best match to the data, and thus which 
one would be used for subsequent analysis, we constructed 
a xt estimate for each simulation by comparing a histogram 
of the residuals in R — i' from the polynomial fit to the data, 
with a histogram of the residuals from a polynomial fit to 
the simulation. The histograms were all constructed in an 
identical manner: we placed the residuals in bins of width 
0.05 mags starting at -0.5 and ending at -1-0.5. Because we 
used a random number generator to select the objects that 
make up the binary sequence, we ran each simulation 1000 
times to obtain a mean value for xt i thus reducing the noise 



for the determination of the most likely binary fraction. Fig- 
ure 13 shows a plot of mean xt against binary fraction. 
Clearly there is no sharp minimum in the value of xt, but a 
broad minimum is centered on 50%. We present one realisa- 
tion of the simulation for a binary fraction of 50% in Figure 

m 

Comparison of panels (a) and (c) in Figure 0] indicates 
that the combination of variability and binarity are not able 
to account for the spread in the PMS. In Figure |5] we show 
histograms of the residuals about 3'^'^ order polynomial fits 
to the observed PMS and the simulated PMS (solid and 
dotted line) for the 15 < i' < 18 region of the CMD. Again, 
it is clear that the combination of binarity and variability 
on timescales of ~ 4 years, are not sufficient to explain the 
spread in C-M space for this PMS. With such a poor match 
to the observed spread it is clear why Figure |3 displays no 
sharp minimum in xi , but rather a broad minimum centered 
on 50%. A binary fraction of 50% maximises the size of the 
spread, for a given estimate of variability. The assumption 
of a 5 = 1 binary population also increases the value of xt 
as it introduces a double peak to the spread which is not 
present in the data. 

The FWHM of the spread in residuals in (R — i') about 
the polynomial fit to the observed PMS, shown in Figure |S]is 
approximately 0.3 mags. This is consistent w ith the width 
of the PMS observed by ISherrv et al. f2004l) in V - I for 
more massive members of the same group. 



3.2 Cep OB3b 

In this case we do not have the benefit of such a large spec- 
troscopic sample as in the a Ori group. As a result, we have 
no estimate of the likely contamination from field stars, or 
the true extent of the PMS region in C-M space. There 
have, however, been a number of studies that have identi- 
fied likely low-mas s members using a variety of techniques. 
lOeura et al. r2002^ used Ha spectroscopy to identify classi- 
cal T-Tauri stars (CTTSs) in the vicinity of bright rimmed 
clouds, and found 33 likely members of Cep OBSb, of which 
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Figure 4. CMDs for the a Ori young group showing: (a) the PMS selected objects as asterisks; (b) the fitted single star and binary 
sequences for a binary fraction of 50%; (c) the simulated PMS for the same binary fraction. 
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Figure 5. Histograms of the residuals in R~i' about polynomial 
fits to the observed PMS (solid line) and the simulated PMS 
(dotted line) for the a Ori PMS sample (binary fraction 50%). 

16 are identified in our catalogue. iNavlor fc Fabian (199^1 ) 
used ROSAT observations to identify 56 X-ray sources to- 
ward Cep 0B3b using both HRI and PSPC observations. 
We have cross correlated their X-ray catalogue with our op- 
tical catalogue, matching the brightest star within a radius 
of 14" for the PSPC positions and 7" for those from HRI. 
We find 21 PSPC objects correlate with objects in our cat- 
alogue, and 14 HRI objects. We reject all objects that lie 
in the Galactic background region on the left of the CMD 
from further use in this st udy as they are likely non-members 
iBurningham et al. 2005)) . identified by chance correlations 
with our catalogue. We accept as likely members those X-ray 
sources that correlate with objects in th e expected PMS re - 
gion of the Vi' CMD shown in Figure ISl lPozzo et al. (20031) 
used radial velocities and Li I absorption to identify both 
CTTSs and Weak T-Tauri stars (WTTSs) in Cep OBSb, 
and found 5 CTTS members and 5 WTTS members. Of 



these, all of the CTTS s and all-bar-one o f the WTTS are 
found in our catalogue. iPozzo et al. (2003^ were also able to 
rule out membership for a number of objects and these are 
also indicated on the CMD. Figure |S|shows the CMD for the 
first epoch data, with membe rs from each survey overlaid, 
along with non-members from lPozzo et al~(2 003). Since we 
are unable to make a reliable PMS selection beyond those 
members selected by previous authors, we simulate the pho- 
tometric spread for these 49 likely members. Since more than 
half our sample has been selected on the basis of X-ray activ- 
ity, it might be argued that our sample is biased in favour of 
objects with the greatest rotational variability: WTTSs. On 
the other hand, WTTS variability is smaller in magnitude 
than non-periodic CTTS variability. As such, our simulation 
should still provide an indication of the contribution from 
variability to the observed spreads. 

We have verified the variability of our member sample 
for Cep 0B3b in the same manner as for the a Ori sample. 
Figure |7|shows the uncertainty in the {V — i') differences as 
a function of V (1st epoch). It is clear that the uncertainties 
start to increase dramatically fainter than about V = 18, 
and step up slightly at brighter than V = 13 for reasons 
similar to those described in Section 13.11 As such we have 
calculated the RMSs for objects in the range of magnitudes 
14 < V < 17.5. The RMS of the diflerences between the two 
sets of observations indicates that the likely members are 
significantly more variable than the rest of the sample. The 
RMS of the differences in V are 0.07 for the total sample, 
compared with 0.11 for the likely members. In (V — i') the 
RMSs are smaller: 0.03 for the total sample; 0.07 for the 
likely members. A similar result is obtained if we attempt 
to make a photometric PMS selection. This indicates that 
a high proportion of objects in this region of the CMD may 
also be PMS members of Cep OB3b. 

The member sequence was simulated in the same man- 
ner as the PMS sample in the previous section. As before, 
we were able to neglect the effect of diffe rential reddening in 
our simulation, but for a different reason. IPozzo et al. (2003^ 
found that the reddening vector in a V/V — i' CMD lies 
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Figure 7. A plot of the uncertainty in A(V' — i') against V for 
the FoV in Cep OB3b. 
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nearly parallel to the PMS for sight-lines toward Cep OB3b. 
As such the differential reddening is unlikely to add any 
spread to the observed sequence. The only difference be- 
tween the two methods is that we have used a different bin 
size and range when constructing the histograms used for 
determining Xv- Because the sample size is smaller we used 
a larger bin width (0.1), whilst extending the range of the 
bins (-2.0 - +2.0) to include some larger residuals. As can be 
seen in Figure |H| a binary fraction of 25% gives the closest 
match to the data. 

The fitted member sequence for this binary fraction is 
shown overlaid on CMD (b) in Figure |2l whilst the simu- 
lated member sequence is shown on CMD (c) of the same 
figure. In Figure 1101 we have plotted a histogram of the 
residuals about the fit to the data and the simulation for 
the U < V < 17.5 region of the member sequence for one 
realisation of the simulation, using a binary fraction of 25%. 
Inspection of Figures Inland 1101 indicates that, as was seen 
in the previous section for the a Ori young group, the com- 
bination of binarity and variability on a timescale of ~ 1 
year is not able to explain the spread of members in C-M 
space. 



Figure 8. The mean for the simulated member sequence com- 
pared to the observed member sequence plotted against binary 
fraction, shown as a percentage for likely members of Cep OB3b. 



4 DISCUSSION 

We have shown that photometric variability on timescales 
of 1-4 years is not able to account for the spread in C-M 
space occupied by either of our sequences. If we assume 
that the influence of variability and binarity are added to 
the underlying distribution in quadrature, we can estimate 
the proportion of the observed spread that is accounted for 
here and that which remains unaccounted for. In the case 
of the a Ori young group it appears that short term vari- 
ability can account for about half of the observed spread in 
Ri' C-M space. The RMS of the residuals about the poly- 
nomial fit is 0.084 mags for the observed PMS, and 0.057 
mags for the simulated sequence, which leaves 0.062 mags 
unaccounted for. Short term variability can only account for 
about 1/20 of the spread in Cep OB3b, where the RMSs of 
the residuals about the polynomial fits are 0.20 mags for 
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Figure 9. A CMD showing the fitted single star and binary member sequences as lines of red asterisks for a binary fraction of 0.25 for 
the Cep OB3b. 
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Figure 10. Histograms of the residuals in V ~ i' about poly- 
nomial fits to the observed member sequence (solid line) and the 
simulated member sequence with a binary fraction of 25% (dotted 
line) for Cep OB3b. 
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Figure 11. A plot of Ai' against A(i? - i') for the 15 < i' < 18 
region of the PMS in a Ori. The points are plotted with error 
bars. 



the observed member sequence, and 0.075 for the simulated 
sequence, leaving 0.185 mags unaccounted for. 



4.1 The nature of the variability 

Figures llll and ll2l show the measured differences in colour 
and magnitude for the objects in our final sample for each 
of our regions of interest. The distribution of C-M shifts is 
clearly different for the two regions, which is not surpris- 
ing considering the different colours used for each case. The 
diagonal distribution of the Cep 0B3b differences demon- 
strates that the variability seen here will only tend to move 
objects up and down an isochrone, and will have little im- 
pact on the apparent age spread. These differences, which 
indicate that objects get bluer when brighter, are consis- 



tent with variability arising from hot or cool spots on the 
surfaces of PMS stars rotating in and out of view. The dif- 
ferences shown in the plot for the a Ori subgroup display an 
uncorrelated distribution, which will tend to displace objects 
across isochrones more. As a result, the simulations carried 
out for this region display a greater apparent age spread. 
Since the o Ori observations are in Ri' , they are less sensi- 
tive to colour changes associated with rota tional variability 
than Vi' observations teerbst et al. 1994) . and so the dif- 
ferent distributions of differences do not necessarily imply a 
different origin. 

Our observations rule out variability on timescales of 
1-4 years as the entire cause for the observed C-M spreads. 
The mechanisms this therefore excludes are rotation of hot 
or cool spots, and short timescale non-periodic T-Tauri vari- 
ability, such as that caused by accretion noise or chromo- 
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Figure 12. A plot of AV against A{V-i') for tlie 14 < V < 17.5 
region of tlie member sequence of Cep OB3b. The points are 
plotted with error bars. 



spheric flaring. We can also rule out variability associated 
with rotation of structures in the disc out to a radius of 1 
A.U as a source for the C-M spread, since the rotation of 
such structures should occur within a year. It is also un- 
likely that rotation of bright structures beyond 1 A.U would 
be responsible for the C-M spread as the temperature in a 
thin disc drops rapidly with distance, such that at a radius 
of 1 A.U the temperature is ~ 100 K. As such the outer disc 
is only likely to contribute appreciable flux in the far-IR. 

The RMS of residuals about the polynomial flt to 
the simulated PMS is larger in Cep OB3 than the a Ori 
young group. Since accreting objects are known to be more 
variable that non-accreting objects, this is consistent with 
the observed incide nce of accretors in the two regions. 
iKenvon et al. found the fraction of low-mass accre- 

tors to be 10 ±5% in the a Ori yo ung group, whilst th e small 
number of members identified bv lPozzo et al. ('200,'^ 1 in Cep 
OB3b suggest about 50% of objects there are accretors. 

Since longer timescale variability may be the origin of 
the remaining unexplained spread we cannot confirm the 
reality of the apparent age spreads. However, what is clear 
is that the actual size of either age spread is smaller than 
that observed, and possibly zero. Returning to the RMSs of 
the residuals about the polynomial fits, we recall that the 
unaccounted for spread has an RMS of residuals about the 
fit of 0.062 mags in the case of the a Ori subgroup, and 0.185 
mags in the case of Cep OB3b. 

It is tempting read much into the observation that the 
larger remaining spread is seen in the association whose na- 
tal molecular would have had the larger crossing time. How- 
ever, it should be noted that in the case of the a Ori sub- 
group the apparent spread, and thus the underlying one also, 
may be larger since we were conservative in our PMS selec- 
tion. As can be seen in Figure 0we have actually excluded 
some members in order to avoid risking significant contam- 
ination by non-members in our photometric selection. The 
excluded members have the same mean variability as the 



members that were selected, so the simulated spread has 
not been altered significantly by their exclusion. 

Although we are confident that our selection 
is not subject to significant contam ination (see 
iBurningham et al. 20051: iKenvon et al. 2005^ ■ there will, 
none-the-less, be some non-members included in the sample. 
These objects will not increase the size of the observed 
spread since our confirmed members span the entire selected 
region of the CMD (see Figure 0. However, it is likely that 
they will reduce the RMS of residuals about the polynomial 
fit to the simulated spread, as any non-members can be 
expected to be less variable than the members. If this is the 
case, then it may still be that the remaining unexplained 
spread in the a Ori is very small. 

It is still possible that the C-M spreads arise as a re- 
sult of some other kind of photometric variability. Such vari- 
ability could arise from variati ons in accretion fiow caused 
by stellar magnetic cycles (e.g. lArmitaee 19951) . with time- 
scales of a few years to decades. Alternatively, much longer 
time-scale variability such as that resulting from accretion- 
driven age spreads (ADAS) could be to blame. Fundamen- 
tally, the infiuence of such long term accretion processes 
on the presence of apparent age spreads would be best de- 
termined through a study of accretion rates across a large 
number of objects in associations that display a range of 
apparent spreads, rather than a longer baseline variability 
study. 



4.2 Apparent age spread and absolute age 

As was briefly discussed out in Section 1.3, the size of age 
spread inferred from a given photometric spread depends on 
the median age of the star forming region in question. This 
should be borne in mind when interpreting any unexplained 
spread in terms of an age spread. For example, consider the 
a Ori young group. If we accept an ag e of 5 Myrs and a dis- 
tance of 350 pc, as we did inlBurinngh am et al. (20*051) . then 
the remaining R—i' spread of 0.062 mags represents a spread 
of approximately 4 Myrs (3.5 - 7 Myrs), based on the sep- 
aration of NextGen isochrones JChabrier fc Baraffe 199?! : 
iBaraffe et al. 2002h at i' = 16.5. If on the other hand we 
adopt the age found bv lSherrv et al. (2004) of 2.5 Myrs and 
the distance they used for that estimate, 440 pc, we find the 
remaining {R — i') spread corresponds to an age spread of 
approximately 2 Myrs (2-4 Myrs). 

In the case of the Cep OB3b association, 
IPozzo et al. (20 03) found the ages of PMS objects to 
range from < 1 Myr to nearly 10 Myrs using isochrones laid 
onto a. V/V — I CMD. As has been shown here, variability 
and binarity can only account for a small fraction of this 
spread. 



5 CONCLUSIONS 

We have used 2 epoch, 2 colour photometry to investigate 
the infiuence of photometric variability on the apparent age 
spreads in CMDs. We have found that the combination of 
binarity and variability on timescales of ~ years cannot ac- 
count for the observed spread in C-M space. We argue that 
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the remaining unexplained spread must either reflect a gen- 
uine spread of ages, or longer timescale variability associated 
with the changes in the accretion flow onto the PMS objects. 
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